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Abstract. The previously derived exact evolution equations for density matrix of 
electron (quantum particle) in phonon field (boson thermostat) are qualitatively 
analysed. Their statistical interpretation is explained in detail, and their main 
symmetry and spectral properties are expounded. In application to the electron's 
random walk, it is shown that these properties certaimly forbid conventionally assumed 
Gaussian long-range asymptotic of the walk statistics. Instead, the exact equations 
imply super-linear dependence of fourth-order cumulant of total electron's path on 
observation time, which signifies existence of 1/f-type low-frequency fluctuations in 
electro's duffusivity and mobility. Physical meaning of this result is discussed, along 
with general origin of 1/f- noise in classical and quantum Hamiltonian many-particle 
systems. 
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1. Introduction 

This paper is direct continuation (second part) of the work [1], where exact equation for 
a quantum particle ("electron") in interaction with harmonic boson thermostat ("ideal 
phonon gas") were reduced to a visual shortened form lightening their consideration, 
and, besides, their relation to 1/f noise (1/f- type fluctuations of electron's mobihty) 
was discussed. 

Here, we start systematic investigation of these equations. First, we will consider 
more carefully their statistical meaning, as it looks from the viewpoint of probability 
theory. Then we will formulate their main principal symmetry properties and spectral 
properties and show that they are incompatible with such "white-noise" statistics of 
the electron's random walk what follows from standard kinetic approximations. The 
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alternative predicted by the exact equations is such statistics which combines white noise 
in velocity of electron and 1/f noise (or "flicker" noise) in its diffusivity and mobility. 

We will finish the paper by discussion of this our concrete result in the framework 
of one more attempt to make our understanding of universal nature of 1/f noise being 
clear also for scientific community. 

2. The model and the problem 

First, let us recollect main results of [1]. We have considered the simplest Hamiltonian 
of quantum particle ( "electron" ) in quantum boson ( "phonon" ) field, 

H = He + Hph + Hint : 

He — ^ , Hph — ^ fiUJk OfeOfe , 

k 

Hint = ^J2^cle^'^ak + c,e-^'^al] , (1) 

and under the thermodynamical limit, when the system's volume O ^ oo , derived the 
following exact hierarchy of evolution equations: 

A„ = - VVx An - i J2 {ook - kV) A„ - 

k<=K„ 

- (thV^o)-' Yl ^k C t^^'"'^' - ^"'""^'(1 + ^k') ] + (2) 

Here 

Aq = Ao(t,X, F) is density matrix of the electron, in the coordinate 
representation, and A„ = A„(t, X, F, i^r„, E„) are definite functions characterizing 
statistical correlations of electron with quasi-classic (coherent-state) ampfitudes 
( = Ak = Qk exp (ikX) , , ^4^ = Al = al exp {—ikX) ) of n phonon modes marked 
by Kn — {ki . . . kn} and and E„ = {ai . . . an} , with aj — ± (or ±1 ); 

Kn^i = Kn Q k and E„_i = E„ in the third term on the right, and 
Kn+i = Kn © q and E„+i = E„ © a in the last tern ("collision integral"); 

X is electron's coordinate, and Y is spatial variable conjugated with electron's 
momentum p , so that V = —i{h/m)VY is electron's velocity operator, and the Fourier 
transform 

j exp {-ipY/h) An d?Y/{2T:hf 

leads to the Wigner representation; 

A^fe = [ex^{tvjJk/T) — are equilibrium occupancies of phonon modes; 

~ '^fc' ^k ~ '^k , 

and flo is some fixed (formally arbitrary) volume. 
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Initial conditions to Eqs.2 can be taken in the form 

An{t^O)^SnoK (3) 

(one particular reasonable choice of Aq = Aq{X,Y) was mentioned in [1]), or 

A^{t = 0) = Wo{X) A^^(y, K^, E„) (4) 

with being stationary (equilibrium and spatially homogeneous) solution of Eqs.2 
(normahzed to Ao^(r = 0) = 1 ). 
Foe more details, please, see [1]. 

Notice that derivation of Eqs.2 there does not touch the purely electron's part of 
full quantum Liouvile (von Neumann) operator, therefore the resulting equations can 
be trivially generalized to electron under action of an external force / (switched on at 
t = 0) hj mere adding terms {—fY/ih) A„ to right-hand sides of (2). 

Our main problem is qualitatively complete and quantitatively correct analysis of 
large-time evolution of electron's position (and total path) probability distribution, 

w{t,x) - Ao(t,x,r = 0) , 

if initially the electron was located in some finite space region, e.g. in vicinity of X — Q . 



3. Statistical contents of the evolution equations 

To perceive more carefully statistical meaning of the functions A„ ( n > ) , let us 
return to their basic definition, - i.e. formula (21) from [1], - and rewrite it in equivalent 
form 

F-\A)P{A)dA , (5) 

where P{A) = P{t, X,Y, A) is the quasi-probability density distribution of all the 
phonon amplitudes A = {A^, A'^} , f .. .dA is integral over all them, and 

The equivalence directly follows from the expansion (20) in [1] (along with relations 
between functions P, Q, Pn and Qn pointed out there). The integration by parts 
turns (5) into 

k keK, 

If all wave vectors of the set Kn are different one from another then this expression 
reduces to formula (38) from [1] . We now will consider general case allowing coincidence 
of some of the wave vectors even in the thermodynamic limit, that is after transition 
from discrete to continuous phonon spectrum. 



lim 



/ 



F(A) 



n 

.k e K„. 



7^0 



dA^" 



^0 



d 



dAl^ 



F{A) } dA (6) 
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Symbolically, this transition reduces to mere re-scaling of the phonon amplitudes, 
so that 



F{A) ^ exp 



(27r)3 ' 



and similar redefinition of the parent creation and annihilation operators, so that their 
commutator turns to 



[afc,aj] ^ {27rf5{k-q) 
(see e.g. [1]) while phonon-related part the Hamiltonian (1) to 

Hph = J hjJka\akdk , Hint = j [cle^''''ak + Ck e''-'^'' al] dk 

with dk = d^k/{2Ti)^ . Consequently, instead of (6) we can write 



{2nf S 



^0 



F{A) \ dA 



(8) 



(9) 



where now J -P(^) {• • •} dA means functional integration. 

The Eqs.6 and 9, together with the mentioned formula (38) from [1], prompt to 
redefine functions A„ as follows: 

A„(i,x,y,x„,E„) ^ A„(t,x,y,x„,-E„) TT —1= 

Then our basic Eqs.2 take form 

An = - y Vx A„ + i ^ (Jfe {Ujk - An + 



keKn 



(10) 



keKn 



d^q 



(27r)3 ' 

where, as above, A„ = A„(t, X, Y, K^, S„) , A„_i = A„_i(i, X, Y, e /c, e (7^) , 
A„_|_i — An+i{t, X, Y, Kn © g. En © cr) . We may write them in compact form 

A = -VVxA + A A , (11) 

with A = {Aq, Ai, . . .} and operator A = A(Y, Vy) unifying all the phonon-related 
operators from right-hand sides of (10). Then Eq.9 yields 

5 



An = /■ P{A) F-\A) n 

I keK„ 



F{A) } dA (12) 
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In the case of non-coinciding wave vectors this reduces to 



or, equivalently, 



A„ = ( diX{t)-X)exp{iYp{t)/h) H A^'it) 



keKn 



in the coordinate representation, or 



An = ( S{X{t) - X) 5{p{t) -p) H A- (t) \ 



(13) 



(14) 



(15) 



in the Wigner representation, where X(t) , p(t) and A^{t) mean electron's coordinate 
and momentum and phonon amphtudes considered as random processes. 

If some of wave vectors from Kn are equal and, besides, their counterparts Ufe = ± 
have opposite signs, then - in contrast to (13), - the braces in (6), (9) and (12) represent 
definite Hermite polynomials of phonon amplitudes. Statistical meaning of the resulting 
expressions becomes clear if one unifies all the redefined functions A„ 's, - with arbitrary 
Kn 's, - into generating functional 



oo ^ „ n 

E /Ann^<^.(^.)^'^. 

n=l o-i...iT„ j = l 



(16) 



Application of Eq.l2 gives 



A = exp 



{2T:fNkZ+{k)z-{k) d^k 



X 



X / exp 



^ j z,{k)Ald'k 



P{A) dA 



(17) 



First exponential on the right here is nothing but inverse of equilibrium characteristic 
functional of free phonon field (i.e. in absence of the electron: 



exp 



{2T:fNk z+{k)z-{k) d^k 



exp 



= ( exp 



I z,{k)Ald'k 
^ 1 z,{k)AlSk 



Pfree{A) dA = 

— free{z} i 



(18) 



free 



where 



F{A) dA 



PfreeiA) = F{A) 

is continuous limit of the discrete measure 

llfk = l[{27rNk)-'exp{-\Ak\yNk) 
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(see [1]). Hence, in the frames of coordinate representation Eq.l7 can be written as 

J'{t,X,Y,z} 



A = A{t,X,Y,z} = 

J'{t,X,Y,z} = 

= ^ S{X{t) - X) exp 

And in the Wigner representation 
A = A{t,X,p,z} = 



J~ free{z} 

lYpit) 



(19) 
(20) 



h 



exp 



/ z,{k)Al{t)d^k 



J='{t,X,p,z} 

•F free{zy 



F{t,X,p,z} = 

= ( 5iX{t)-X)Sipit)-p)exp 



(21) 
(22) 



J2 J ^a{k)Al{t) d'k 



Formulae (17)-(21) give complete and transparent statistical interpretation of generating 
functional of A„ 's and thus A„ 's themselves. It remains to notice that in (symbolic) 
terms of the full density matrix of the system, p (see formula (4) in [1]), 



F{t,X,Y,z} = IVp^exp 



X (r| p \r') exp 



ttk e'''^ d^k 



X 



e-^'^'"' d'^k 



(23) 



where r — X + Y/2, r' — X — Y/2, \r) and \r') are eigenstates of the electron's 
coordinate operator, and Tfp^ is trace over all phonon states. 

4. Generating evolution equation 

The Eqs.lO are useful equivalent of Eqs.37 in [1] for functions Dn representing the same 
correlations as A„ in Eqs.l4 and 15 but in dimensionlcss relative units (obviously, 

A„(X„,E„) = L>„(X„,-E„) II 

keK„ 



after our above redefinition of A„). Corresponding equivalent of the functional 
evolution equation (40) from [1] is 



A = - VVx A + 



d^k Za{k)Lka 



5za{k) 



A + 



(24) 



+ 



d'^k z^{k)Bka A + 



6 



d^k Ak, 



5za{k) 



VVxA + C{z,-) A , 



where now 



Lk,j = icr{ujk - kV) , 
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B,. ^ ^ W'^^'N, - e-'^^'I^N, + 1) ] , (25) 
in 

Correspondingly, the functional J-" undergoes equation 

-VVxJ" + cAz,^\ F , (26) 



dz 



where 



■^0 — J' free ^ J' free ~ 



Y, jSk.Ak)L,,^ + (27) 

Sigma 



(k) 



with 



BL ^ e-^^/^ (28) 

Notice that the evolution operator Co does not include the phonon mode occupancies 
(which instead must appear in initial conditions for J"). 

5. Symmetry properties of correlation functions and their evolutions 

According to previous section, both the hierarchy A = {Aq, Ai, . . .} as the whole and 
the functional (16) are equivalent to full density matrix of our system. The hermicity 
of this density matrix means that 

A4-r,-E„) = A;(r,E„) , (29) 

in the coordinate representation, or 

7^A„(y,E„) = A;(-y,-E„) = A„(y,E„) 

Here and below we write out only arguments what are under current attention. In terms 
of the generating functional, 

nA{Y,z,} = A*{-Y,z_,} = A{Y,z,} , (30) 

where * means complex conjugation of coefficients of the functional (but not its 
argument Za{k) ). Indeed, one can verify that evolution operator in Eqs.lO and 24 
satisfies 

n[-vvx + jC]n-^ ^ [-vvx + c] (31) 

Hence, if the hermicity property is satisfied at any one time moment (e.g. when initial 
conditions (3) take place) then it keeps at all other time moments. In the Wigner 
representation the same property looks as 

?^A„(p,E„) = A;(p,-E„) = A„(p,E„) , (32) 
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with similar changes in (30) and (31) and with V =^ p/m . 

Further, if the system's Hamiltonian expressed by (1) and (8) is invariant in respect 
to time inversion, - which is the case when 

uj-k = (^k , c-k = cl , (33) 

- then 

e[-vvx + A]e-^ = - [-vvx + a] , 

e [ -vvx + c]e-' = - [ -vvx + c] , (34) 

where © is time reversal operator defined by 

eA„(y,ir„,E„) = A„(-y, -E„) 

eA{Y,z,(k)} = A{-Y,z_,{-k)} (35) 
in the coordinate representation and by 

e A„(p, E„) = A„(-p, -S„) , 

eA{p,z,{k)} = A{-p,z_,{-k)} (36) 

in the Wigner's one (in any time-reversed process, all momenta of particles and 
quanta are inverted while their creations are replaced by annihilations and vice versa). 
Hence, if A„(i, X, Y, X^, E„) and A{t, X,Y, Zcr{k)} are some solutions of Eqs.lO and 
24, respectively, then © A„(to - t, X, Y, K^, E„) = A„(to - t, X, -Y, -K^, -E„) and 
© A{to — t, X,Y, Za{k)} = A{to — t,X, —Y,z-a{—k)} also are their solutions. As the 
consequence, the equilibrium stationary solution of evolution equations satisfies 

©A^« = A:^ (37) 

The properties (31) and (34) are characteristic principal properties of (quantum) 
Liouville operators. This is quite natural: relations (19), (21) and (23) do show that 
the evolution operators —Wx + A and —Wx + ^ are nothing but particular 
representations of the full system's Liouville operator, i.e. [H, . . .]/ih. 

This means that operators —Wx + A and —Wx + ^ have the same 
spectrum as [H,...]/ih do. That is all their eigenvalues are purely imaginary, since 
spectrum of [H,...]/ih consists of {Es — Es')/ih, with Es being eigenvalues of 
the full Hamiltonian: H\S) — Es\S) . According to relation (23), the corresponding 
eigenvectors of —V"Vx + ^ and thus —V"Vx + A are determined by expression 



Ass'{X,Y,z} = ^i,Uz} Tiphexp J z_{k) a^e'''' d'k 



X 



X (r| 15") (S*'! \r') exp 



j z+{k) 4 e-'^""' d^k 



(3J 



Clearly, the functionals introduced by this expression form definite complete 
set of cigcnfunctions of our evolution operators. Then the imaginary character of 
corresponding spectrum means that these operators can not have real eigenvalues or 
complex ones with nonzero real parts. Next, we will consider most principal and 
important consequence of this circumstance, while Appendix A contains a its simple 
visual illustration. 
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6. Joint characteristic function of total electron's path and instant system's 
state, and their cumulants 

Let initially the electron is located in some finite spatial region around the coordinate 
origin, so that, for example, 

A°(X,y = 0) = W^{X) = (27rxg)-=^/^exp(-XV2xg) 

in (3) and (4), respectively, with some reasonable uncertainty xq (e.g. Xq = fi? /AmT 
[1]). In the first case, the electron is suddenly thrown in the phonon gas just at initial 
time moment t — Q and therefore at i = has no correlation with phonons. In the 
second case, at i = it already possesses all equilibrium correlations with phonons 
correspondingly to momentum and energy balance with them. 

This difference, however, must become unimportant after sufficiently large time, 
t ^ Tq , where t ^ Tq is characteristic time of relaxation (thermalization) of electron's 
momentum probability distribution. Therefore at t ^ Tq ,- or, formally, at t ^ oo , - we 
can concentrate on statistics of electron's random walk, i.e. on probability distribution, 
W{t,X) , of its coordinate X — X{t) . At that, since the latter unboundedly spreads 
in space, - so that {X'^ (t)) / x^, ^ oo , - we can treat X{t) as total electron's path 
accumulated during time interval (0, t) . 

Then, as usually, it is convenient to investigate firstly characteristic function of the 
path distribution, by considering Fourier transforms 

A„(t,i«;,y,i^n,En) = /exp(ifi:X) An(t,X,y,i^n,E„)d3x , 
A{t, in, Y,z} ^ f exp {inX) A{t, X, Y, z} d^X , 
so that the evolution equations take form 

A„ = iKVA^ + Lka,A„ + 
keKn 

+ J2^ka,K-i+ Yl [AgaK+id^q (39) 

(where, of course, K^-i = Q k , E„_i = E„ (7^ , K^+i ^ ® q and, 
E„+i = E„ © (7 ), - or shortly 

A = [inV + A] A , (40) 

- and 

A = iKV A + CA , (41) 

with operators defined in (24)- (25). 
According to Eq.l9, 



where 



F{t, iK, Y,z} = J exp {inX) J^{t, X, Y, z) d^X = (43) 
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= ^exp iKX{t) + iiV{t) Za{k)Al{t) d^k 

with V{t) —p{t)/m being electron's velocity, and 
i = mY/h 

Thus, F{t, IK, y, z} in fact is full characteristic function of our system (in the sense 
of the probability theory), while A{t,iK,Y, z} differs from it only by exclusion of 
undisturbed part of self-correlations of phonon modes. At the same time, since 

X{t) = X{0) +f^ViT)dT 

and \X{0) \ ^ xq , function A{t, in, Y, z} serves as characteristic functional of electron's 
velocity on the interval (0, t) . 

As usually, it is reasonable to consider this function in terms of various cumulants 
and their generating functional: 

C{t,iK,^,z} = In A{t,iK,Y,z} = (44) 



E 



j\Un\ 



/ ^o-i(^l)-- - / ^<Jn{kn) Cjin{t, Knj'^n) ) 



j+l+n > 

C,,„(t,X„,EO = Ux^{t)V\t) n Allit) 



s = l 



where J^. . . = J ■ ■ ■ d^k , Kn = {ki . . . A;„} , S„ = {ai . . . o"„} , the double angle brackets 
denote (j + l + n) -order cumulant, i.e. purely irreducible part of (j + l+n) -order mutual 
correlation, of the enveloped multipliers, and equity of some of indices j,l,n to zero 
means absence of corresponding multipliers (only in special case, when j = / = 0, n = 2 , 
Coo2 is difference between full cumulant and its "seed" value determined by In J^jree ) • 
In particular. 



C{t,iK,i,0} = In Ao(t,iK,Y) 
(in)^(tO' 



(45) 



j+l > 



j\l\ 



{{X^(t)V^(t))) 



describes most interesting for us full joint statistics of electron's velocity and total path 
in themselves. 

Notice that in terms of cumulants Eq.41 becomes nonlinear: 

dC 



+ ^ ^ d^k Za(k)Lk^ 



d k Za{k) iak 



6za(k) 



VC + 



(46) 



sc 



+ 



6za(k) 

+ E/^^^ ^MBka + ^jd'^k Ak, 
with quadratic nonliearity. 



(k) ' 
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7. Long-time asymptotic of the cumulants and electron's 
diffudivity /mobility low-frequency fluctuations 

Now, let us consider long-range time behavior of the electron's path and related 
cumulants at i/ro — >■ oo . 

Our sole assumption will be that electron's interaction with the phonon thermostat 
ensures fast enough, - i.e. integrable, - relaxation of the electron's velocity correlation 
function (( V{t)V{t') )) and thus diffusive character of the random walk, in the sense 
that 

C2oo(t) = {{X\t))) ^ 2Dt + const (47) 

Undoubtedly, this is the case under suitable choice of electron-phonon couplings |cjfcp 
and phonon frequencies ouk ■ In particular, without essential loss of generality we may 
take 

Ck = c{\k\) , ujk = uj{\k\) , (48) 

with c(|A;|) = c*(|A;|) > (phases of Ck 's anyway are of no importance). This guarantees 
spherical symmetry of the random walk. 

A suitable choice can be recognized from the simplest approximation of exact 
equations, i.e. the "kinetic equation", - see Eq.58 in Sec. 5. 4 in [1]. In the representation 
under consideration it can be written as 

Ao{t,iK,Y) = [ikV + K]Ao{t,iK,Y) , (49) 
K . . . = - ^ / Aka [Lka - ■ ■ ■ , 

c 

with e = +0 t . 

Provided the mean square of electron's path obeys the diffusive law (47), let 
us discuss possible long-term asymptotic of higher-order path cumulants and other 
cumulants. 

If we were basing on habitual intuition, in turn based on standard kinetic theory, 
we would predict that all the cumulants from (44) and (45) tend to finite hmits, except 
self-cumulants of the electron's path: 

CsirSS) const if l + n> , (50) 
Csooit) = {{X'{t))) ^ sD,t + const , (51) 

where 

^'^llt^^^^^^^^^ ^ {{nt)X'-\t))) ^ const 

t 

Notice that Lfco- is imaginary-valued operator. Addition of infinitely small real positive e = -|-0 to 

it is necessary for artificial breaking of the time reversal symmetry. Otherwise, at e = , as one can 
easy verify, the integral in (49), i.e. K , would turn to zero (see also Appendix A below). 
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and D2 = D (of course, Z^^ = at odd s because of the spherical symmetry). Such 
behavior means that, at any n > and any i]j{t) being either V{t) or some of A1{t) , 

{{m{t)...Vs{t)Vit,)...V{Q)) ^ (52) 
when some of t — tj ^ 00 , 

where the zero hmit is achieved in so fast way, that these cumulants are integrable over 
all tj simultaneously. Orally, all irreducible correlations between present state of the 
system, - characterized by electron's velocity and phonon amplitudes, - and past values 
of velocity disappear with time. Then we can write 

A{t,iK,Y,z} = e^^'""^^ A'{t,iK,Y,z} e^^**^)* A' {00, m,Y, z} , (53) 

with 

CO, m, 0,0} ^ , / n 

XiiK) = lim = > — tttI-k ) < 

^ ^ t->oo t ^ (2s- 1)! ^ ^ 

Correspondingly, the path probability distribution is asymptotically Gaussian: 

W{t,X) -> {47rDt)-^/'^exp{-Xy4Dt) (54) 

Just such type of asymptotical behavior follows from the kinetic equation. At that, 
A(m) is real negative quantity determined by the eigenvalue problem 

X{iK)f{m,Y) = [mV + K]f{iK,Y) 

But, unfortunately, from the viewpoint of exact equations, (39) or (40) or (41), so 
nice behavior is impossible! 

Indeed, if the assumed asymptotic (50), (51), (53) was true, then, evidently, the 
Eq.41 would result in equality 

[iK,V + C] A' {00, in,Y,z} = \{iK) A' {00, iK,Y,z} , (55) 

thus stating that operator inV + C has a real eigenvalue. Then this would imply, - as 
it is easy to see, - that operator —Wx + ^ (and, equivalently, —Wx + ^ ) has the 
same real eigenvahic, 

X{tK) f{X, Y,z} = [ -VVx + C] f{X, Y, z} , (56) 

with corresponding eigenvector expressed by 

f{X,Y,z} = e-'"^ A' {00, iK,Y,z} (57) 

However, this is impossible, since, - as we have shown above, - the whole spectrum of 
the exact evolution operator is purely imaginary. It is useful to add that this reasoning 
stays valid even if the constants in (47) and (51) are replaced by arbitrary sub-linear 
time functions. 

Hence, the assumptions (50) and (51) (at s > 2), and thus (53) and the 
hypothetical asymptotic (54), in spite of their seeming plausibility, all are qualitatively 
wrong. 
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One more proof of incompatibility of these assumptions with exact evolution 
equations comes from the property (34). The latter implies, firstly, that if (57) is 
eigenfunction of the operator —Wx + C, , with eigenvalue A(m) , then this operator 
has also eigenvalue — A(i«;) , corresponding to eigenfunction e~"^^ © A'{oo, ire, Y, . 
Therefore, secondly, starting at i = from initial conditions (e.g. (3) or (4)) what 
satisfy time- reversal symmetry, like (37), we must come to asymptotic which, in contrast 
to (53), with equal rights includes two exponentials, exp [A(i«;)i] and exp [— A(iK)i] , 
and thus is certainly physically senseless! 

Consequently, we have to conclude that in reality some of the ccumulants in (52) 
arc not fast enough decaying, i.e. integral of some of them over all tj from to i are 
unboundedly growing with time t . The natural first candidate for such role is § 

that is 

{{V{t)X\t))) ^ oo (58) 
This implies that 

Thus, the fourth-order cumulant of electron's path grows with time obeying 
a super-linear law. This conclusion is main final result of the present paper. 

We see that long-range behavior of (even equilibrium) electron's random walk can 
not be adequately described with the help of single parameter D , which now determines 
mean square value of path only, but not its higher-order statistical moments. 

Further consequencies of this result, concerning the 1/f -noise, do follow already 
known general relations [3, 4, 6, 7, 8]. Namely, the asymptotic (59) shows that the 
electron's path X{t) behaves as if electron's diffusivity was a random variable, Dii) , 
with mean value D and effective correlation function 

((l)(t)5(0))) = = \ {{V{t)X^{t)V{0))) (60) 

From (58)-(59) one finds that this "diffusivity correlation function" is slowly decaying 
(definitely non-integrable) time function || . 

The cumulant {{V{t)V{ti))) , can not play such role as far as the mean square diffusive law is 
fulfilled, either in the form (47) or {{X^{t))) /t H- 2D = const. 
II 

Formulae like (60) appeared in [3, 4, 6]. They presume that observation of (equilibrium) random 
path starts from nearly equilibrium initial condition like our (4). 

One of possible forms of the diffusivity correlation function was predicted and considered in [3, 4, 6, 8]. 

In principle, generally speaking, this function, being introduced by (60), can be even not decaying but 
instead growing with time. Such possibility was demonstrated and explained from physical point of 
view in [7, 8, 9, 10, 16]. 
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Analogous consideration of our equations in presence of external force (as was 
mentioned above) or, alternatively, application to our results of the generalized 
fluctuation-dissipation relations shows that the diffusivity fluctuations always convert 
themselves into similar fluctuations of (low-field) electron's mobility, with the same 
correlation function (60) % . 

Obviously, the corresponding spectrum of diffusivity /mobility fluctuations is of 1/f 
type"*" - or, more generally, of 1/f type, - showing unbounded growth at /tq 0. 
A concrete form of this spectrum produced by our exact equations will be investigated 
separetely. 

8. Discussion 

As it was realized already in [3, 4, 5, 6, 7, 14], low-frequency 1/f -type fluctuations 
in diffusivity and mobility of charge carriers, - as well as other molecular-size (quasi-) 
particles, - originates from indifference of many-particle systems to a number and relative 
frequency of collisions (elementary acts of interactions, scattering, etc.) of any concrete 
particle. In other words, indifference of many-particle systems to rates of participation 
of a concrete particle in various irreversible and random processes, first of all, in the 
particle's thermalization (e.g. momentum and energy relaxation), diffusion and drift 
under external forces. 

Let, for example, during last 1 second some air particle or charge carrier has 
undergone two times greater (or lesser) number of collisions than it must have "on 
average". Such incident will make no influence upon next particle's motion, since any 
detail of its motion becomes forgotten by the system after characteristic time Tq much 
shorter than 1 second. In other words, even strong fluctuation of number of collisions, 
- and hence of rate of particle's diffusion, - does not cause a "back reaction" of the 
system (which would enforce this particle to make "opposite fluctuation" in number of 
collisions during next 1 second and "compensate" the old incident). 

This reasoning shows principal difference of fluctuations in relative frequency of 
collisions (and connected physical quantities) from, say, particle's energy fluctuations 
(which cause back reaction suppressing energy deviations from its mean value during 
time ~ To ) . 

The 1 second in above reasoning can be replaced by arbitrary time interval greater 
than To , and instead of 1 second we can take 10 seconds, 100 seconds, and so on * . 

1 

See [5, 6, 7, 8, 11, 12] and, for the fluctuation-dissipation relations, references mentioned in [6, 7, 8, 13] 

and in preprints [Yu.E.Kuzovlev, arXiv: 1106.0589-,, 1108.1740]. 

+ 

For details and examples see e.g. [4, 6, 7, 8, 9, 10, 11]. 

Strictly speaking, this interval can be bounded above by some characteristic "non-ergodicity time" , 
which, however, is very large for real, even closed, many-particle systems [10]) and infinitely large for 
open or infinitely many-particle systems. 
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At that, because the system constantly forgets history of the particle's collisions, it, 
figuratively speaking, can not distinguish what part of collisions is "average" and what 
is "deviation from average". Hence, it allows the latter to be as fast growing with 
time as the former. This characteristic law, N{t) — {N{t)) oc {N{t)) oc t (with 
N{t) standing number of collisions) just means that spectrum of fluctuations of instant 
relative frequency of coUisions (and thus of rate of diffusion, etc.) at low frequencies 
/ <S I/tq has nearly 1/f form (differing from it by some logarithmic factors). 

Additional explanations of so simple origin of 1/f -noise can be found in [1, 8, 9, 10, 
11, 12, 13, 15, 16]. 

Quite obvious conclusion from the aforesaid is that theoreticians should honestly 
investigate actual Hamiltonian dynamics of many-particle systems and their dynamical 
chaos, - basing on exact equations of statistical mechanics, - instead of thoughtless 
demonstration of BernouUi's "art of conjectures" . 

In essence, the same was claimed by N.Krylov [17] as far ago as in 1950. One 
of important results of [17] is statement that even in "good" Hamiltonian systems 
(possessing excellent mixing properties) relative frequencies of a given sort of events 
on different phase trajectories, - i.e. in different experiments, - generally differ one 
from another and thus from ensemble average, regardless of duration of experiments. 
Therefore, - as N. Krylov underlined, - it is impossible to describe actual randomness 
of physical systems in terms of any a priori introduced "probabilities of events" (e.g. 
collisions). In other words, "probabilities of events" have no certain values. 

Hence, all probability-theoretical models of noise inevitably lose fundamental (in 
the above explained sense) 1/f -noise. That is why we have to investigate dynamical 
models without any a priori truncation of their exact equations (for any truncation, - 
like one leading to the "kinetic equation" (49), - acts as artificial introduction of some 
a priori probabilities). 

One more important Krylov's result is that in statistical mechanics statistical 
correlations between particles or events do not necessarily say about some real 
"physical" (cause-and-consequence) correlations or connections between them. Our 
above treatment of collisions' number fluctuations well illustrates this fact. Thus it 
would be wrong to interpret the 1/f -noise under discussion as manifestation of some 
physical long-living correlations (e.g. between charge carriers) or some "slow relaxation" . 
There are only purely statistical correlations, and they manifest only so rich randomness 
which can not be adequately imitated by traditional probability-theoretical (kinetic and 
stochastic) models. 

The situation with 1/f -noise in semiconductors [18, 19] serves as remarkable 
example (see [6, 8, 12]). On one hand, many experiments give strong evidences that 
main source of this noise is (independent) 1/f -type fluctuations of mobilities of charge 
carriers resulting from their scattering by lattice vibrations, i.e. phonon gas [18, 19]. 
On the other hand, the traditional way of thinking enforces to switch attention from the 
electron-phonon interaction and scattering itself to the phonon subsystem in itself and 
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suggest 1/f fluctuations of occupancies of phonon modes (see e.g. [19] and references 
therein). Even in spite of clear contradiction between such idea and kinetic theory 
which establishes quite definite (finite) relaxation times of phonon modes! jj . In fact, 
however, low- frequency 1/f -type fiuctuations observed in phonon systems themselves 
(e.g. fiuctuations of internal friction in quartz crystals or intensity of light scattering 
by them) come from fiuctuations in relaxation rates of phonon modes there (see [13] 
and references therein), - in accordance with the above reasonings, - but in no way from 
fiuctuations of phonon numbers. 

In our Hamiltonian model, considered in this paper, perturbation of any particular 
phonon mode by the electron vanishes under the thermodynamic limit. In other words, 
any finite sequence of ^ t/ro electron's scattering events during finite time t involves 
only vanishingly small portion of total number of phonon modes. Therefore, all their 
occupancies do not fiuctuate at all, staying mere random constants. At that, since they 
are independent one on another, all effects of their randomness are suppressed by the 
infiniteness of their density in k -space (any finite volume there includes infinitely many 
modes). Hence, equations of our model describe just what is most interesting for us and 
for applications, i.e. 1/f -noise resulting from uncertainty of rate of electron's scattering 
by phonons. 

Of course, real samples under measurements are finite and rather small, but one 
should remember about very many degrees of freedom in their surroundings which 
inevitably somehow interact with any sample. 

9. Conclusion 

To resume, we have shown that if an interaction of electron (quantum particle) with ideal 
phono gas (harmonic boson thermostat) enforces the electron to move like Brownian 
particle, - so that mean square of its path obeys the diffusion law, {X^{t)) oc t , - then 
such interaction certainly supplies this motion with 1/f (or "flicker") low-frequency 
fluctuations in electron's diffusivity and thus mobility. 

In essence, this our result was deduced from most general properties of exact 
statistical-mechanical equations of the electron-phonon system only. This fact prompts 
that our result can be extended to a wide class of different systems, in the form of 
some general theorem (which surely will include the molecular random walk in fluids 
considered in [7, 10, 11, 16]). This seems very attractive task. But not less intriguing 

tt 

Thus authors of such idea are even more "revolutionary" than I am! I do not suggest that results 
of kinetic theory are invalid, stating only that they are incomplete, and that, to supplement them 
properly, one should return from habitual probability-theoretical scheme 

Stochastic 

Probabilities ^ model Noise 

to canonical statistical-mechanical scheme 

Hamiltonian 

model Noise Probabilities 
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task is concretization of our present result up to exact description of shape (frequency 
dependence) of the 1/f -type spectrum, along with its adequate enough quantitative 
estimates. I hope to contribute to such work. 

Appendix A. Imaginciry spectrum of the evolution operator and falsity of 
the kinetic equation 

Notice that properties (31) and (34) remain formally valid under truncations of the 
evolution equations. Therefore character of spectrum of the evolution operator remains 
unchanged. And, in order to feel it, we can for simplicity make truncation already at 
n — 2 , by formal putting A2 = , as was done in [1] when deriving the "kinetic 
equation", Eq.58 in [1]. Besides, for more simplicity, we will confine ourselves by 
spatially homogeneous eigenf unctions. 

Thus, let us consider the eigenvalue problem 

^^o(r) = -l a/cne--/^-e--/^)A,(r,^.,.)^, 

XA,{Y,k,a) = }^ac,-[N,e^''''/^-{N, + l)e-''''^/']Ao{Y) + 

+ ia{ujk - kV)Ai{Y,k,a) , (A.l) 

and assume that it has solution with a real non-zero A . Prom viewpoint of 
the kinetic equation this seems undoubted, because kinetic equation arose in the 
same approximation, and spectrum of its X -independent eigenfunctions (spectrum of 
momentum and energy relaxation) certainly is real. 

Excluding Ai from this pair of equations and going to the Wigner representation, 
we come to 



J \ck\'^{Sx{uJk- kV +^) X 



AAo(p) = (A.2) 

X [ Nk Ao(p - hk) - {Nk + 1) Ao(p)] + 
+ dxiujk -kV - —) X 

X [(Afc + 1) Ao(p + hk) - Nk Ao(p) ] } "^'^ 



/ 



(27r)3 

[V^A(pb')Ao(p') - V^A(p'b)Ao(p)]dV , 



where 



+ 



27r r 

Vx{p\p') = ^ j \ck\^5{p-pt-hk)[Nk5x{uk + hk^/2m-kV) 
+ (Afe + 1) 5x{uk - hky2m + kV) ] -0^ (A.4) 
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At A — )■ +0 operator on the right in Eq.A.2 turns to the kinetic operator from 
the Eq.58 in [1], that is operator K from the same Eq.49) above, whose spectrum is 
real non-positive. At finite A > , evidently, the operator in Eq.A.2 has the same 
structure as general kinetic operators, therefore it also is non-positively defined. Hence, 
Eq.A.2 can not have non-trivial solution with A > . If A < then signs of both 
sides of A < do change simultaneously and thus remain opposite. Therefore A < 
also can not give the solution. Clearly, addition to A of imaginary component does not 
change the situation. Consequently, all non-zero eigenvalues in the Eqs.A.l are purely 
imaginary. 

This well illustrates that the kinetic equation, in spite of its usefulness, is in principal 
contrast with its parent evolution equations. 
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